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PART It THE CASE OF N0N81N6ULAR THR08T 
Blieser l^lndlar and Frmk Neuaum 

I. Introduction 

The literature on fllghtpath optlodaatlon la extenalve. It can be 
classified according to paths in the vertical plane, the horisontal 
plane, and the three-dimensional space; It can be further classified 
eccordlng to the type of aircraft and mission, and the performance Index. 
We consider minimum-fuel, constant-altitude fllghtpaths of a transport 
airplane in the terminal area. 

Host of the papers on fllghtpath optimization in the horizontal 
plane consider minimum tlme.^'^ In Refs. 1-3, the velocity Is constant; 

In Refs. A-7 the velocity Is a state-variable as In our case but the 
assumptions, constraints, and numerical results correspond to a supersonic 
fighter aiicraft. Thus, although there is contact with our results, 
there Is no overlap. References 8 and 9 consider minimum-fuel, horizontal 
rocket turns, but since the mass Is variable, the results are quite 
different. An overview on fllghtpath optimization Is given in Ref. 5. 

Our objectives in investigating a minimum-fuel landing problem were 
to gain Insight into the characteristics of minimum-fuel fllghtpaths by 
analysis and computation, and to use these results to Improve the on-line, 
fuel-efficient capture algorithm of Ref. 10. Details of the refined 
algprlthm (in the horizontal plane only) are reported in Ref. 11. 

Following statements of the problem and of the necessary conditions 
in Secs. II and III, respectively, the optimality of strai^t-line and 
circular fllghtpath elements (which the algprithms of Refs. 10 and 11 
use) is Investigated in Sec. IV. It is shown that a straight-line 
seg^nt can occur only at the beginning of a minimum-fuel fllghtpath; 
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this also ravaalad tiow aucii patto can ba eoaputad. Cirealar paitha ara, 
la paaMalt aM: Caal>apclay.. Caaputatim of rap a aaa ata tiva axtraaaia 
la diaeuaaad in Sac. V. fha raaultlng axtranal fll^tpatba can ba 
groiqmd in tluraa catagoriaat oaaatinn anttanala . daealarating wil:b 
aero thrust throughout; short«ranaa tumina axtranala » where the initial 
and final positions* but not the headings* are relatively near* say* 

1-3 n. mi.; and long-range extremals * characteriaed by a possible Initial 
turn* followed by a long (say* 5-15 n. mi.) almost straight arc and 
ending with a possible final turn. Since the global optimality of some 
of these extremals was suspect, they were checked against near-optimal 
flightpaths produced by the algorithm of Ref. 11. nils comparison, 
discussed in Sec. VI, established the existence and approximate location 
of Darboux points (beyond which the extremal ceases to be globally 
optimal; see Bef. 12). 


II. Problem Statement 


The point-mass equations of motion in the horizontal plane are 


V cos <|) 

(1) 

V sin Ip 

(2) 

-gu/v 

(3) 

(T - D)/m 

(4) 


Here, x and y are the coordinates of the horizontal plane, ^ the head- 
ing angle measured counterclockwise fr^m the x-axls* v the velocity, 
g the gravitational constant, and m the mass; the control variables 
are the thrust T and u, where u - tan ^ and ^ is the bank-angle* 
positive with right wing down. The drag D is given by 
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D • kjV* + kjd + u*)/v* 


(5) 


Where k^ and k^ are conatanta (i*e*( they are aaeuaed to be independent 
of V for the low velocltlea in the temlnal area) . Theae equatioi» 
were uaed in Refs. 4-7 and are derived by aaauning zero wind velocity • 
constant mass, coordinated turns, and a snail angle of attack which is 
automatically adjusted to maintain horizontal fll^t (see Appendix A). 

Constraints and final states are 


|u(t)| <u„ 


T^n 2 ^ W 


v(t) < V 


max 


( 6 ) 

(7) 

( 8 ) 


x(tf) - y(tf) ■ 0 , i|/(Cf) - k2ir , k - 0, ±1, . . v(tf) - Vf 

(9) 

The cost Integral to be minimized is the fuel consumption 

/ tj 

Icj + CjT(t) + C 2 T^(t)Jdt t C£ ■ const >0, i*0, 1,2 

0 

( 10 ) 

where t£ is free. The terms Cg and C 2 T^ in the fuel-flow-rate are 
often neglected in the literature. For the case considered here, Cg, 
the fuel-flow-rate at zero thrust, is not negligible. The term CgX^ is 
small but significant: when C 2 * 0, the optimal thrust is discontinuous 

and its intermediate values are singular. By changing the units of the 
C£, the cost integral J can be Interpreted as a combination of the 
cost of time and fuel (the operating cost). The time-optimal problem, 

Cp ■ 1, Cj ■ Cj ■ 0, has been treated in Refs. 1-7 and is not considered 
here. 
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In suDBMryi elw problna Is to dctsmlno the controls u(e) snd T(t) 
mi%A Mrs^nifiBdlaa stoto tmioctocv Cron s& sidilrsarv ■»»»^ 

St the tlBM t ■ 0 to ths flnsl ststs of Bq. (9) st s frss tlM t ■ t|« 
subject to Bqs. (l)-(8)i so ss to nintotss tlis fusl <wnsusptlon CBq. (10)]. 

For our gensrsl snslytls» vs maks the following two Msunptions, 
which are easily satisfied for our (and Indeed for aoet) ntnerical 
values. 

Assumption 1 . For the applicable range of velocities, the thrust 
that equals drag is intermediate. This implies that the two velocities 
for which T ^av “ D are outside the applicable range. It also implies 
that is less than the minimal drag with respect to velocity. 

^mln' general case, can be negative, and we make 

Assumption 2. 

Assumption 2 . is such that the fuel flow rate at T > 

is positive and ^mio ^ -c^/(2c 2). The first part eliminates considera- 
tion of gliding flightpaths with shut-off engines; all such paths are 
fuel-optimal and trivially satisfy the necessary conditions embodied in 
the minimum principle. The significance of the second part is discussed 
in the next section. 

Numerical values used in this study correspond to a 150. 000- lb Jet 
transport at sea level at velocities not below 150 knots. 


V “ 250 knots , v, > 180 knots 


T - 30,000 lb . 
max 


min 


0 


A > 30® , tan ^ ■ 0.577 

in n in 
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kx ■ 0.08 Ib/knot^* kx ■ 2.127 x 10* lb knot* 

Cj ■ 0.808 lb O 0 c“* » Cj ■ 1.507 x 10”* ooc”*, Cj ■ 5.4 x 10“** lb“* sec 

To get an Idea on the percent^ contributlona of the terns c^T, 
and C 2 T* In Eq. (10) to the total fml-flow^rate, aaaune a flight along 
a stralglt line at a constant velocity of 250 knots (at T <■ D ■ 8.403 lb). 
The contributions are then 38Z for Cq» 60Z for c^T. and 2Z for C 2 T*. 
Assumption 1 is easily satisfied: we find that Tma» ■ D occurs at 
83 knots and 607 knots, and that is 8,250 lb at 227 knots. 

Evidently, Assumption 2 is also satisfied. 

III. The Necessary Conditions 

We employ the minimum principle. The Hamiltonian is 
H ■ Xq(co + CjT + CjT*) X^v cos i|* + XyV sin i)> - X^gu/v 

+ Xy(T - kjV* - k 2 (l + u*)/v*l/m + n(v - v^air ) (H) 

where n > 0, n(v - = 0 (Ref. 13). The costate variables are given 


by 

Xq ^ 0 , Xo <■ constant (12) 

Xjj « -Hjj ■ ® “* *x " constant (13) 

Xy ■ -Hy ■ 0 •* Xy ■ constsnt (14) 

X^ - -H^ - v(Xjj sin p - Xy cos p) (15) 

\y - -Hy - -Xjj cos p - Xy sin p - X^gu/v* 

+ 2XytkjV - k 2 (l + u*)/v*J/m - n (16) 
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irtMvc - 8R/8x» «tc. Sine* U not pnocHlwd and • 0, no 
Invo tho nMOsaoxy condition 

H s 0 • all t e (0, t^] (17> 


If for Bone t 


l^j»(t) ■ Xy(t) ■ 0 • ■ Xy 


(18) 


then Eq. (17) and Assumption 2 require the vanishing of Xg. Since having 
ell coststes zero Is not optimal . an extremal where Eq. (18) occurs is 
not fuel optimal. Henceforth* we consider only extremals with X, > 0* 
and we normalize the costates by setting Xq ■ 1. 

Minimization of H with respect to T yields the extremal thrust. 


where 


T* 


max 

X 


if 


T > T 
* max 




min 


if 


T < T . 

“ min 


(19) 


X ■ -(Cj + X^/a)/(2c2) 


( 20 ) 


We note that by Assumption 2, Eqs. (19) and (20) show at once that 

Since the minimization of U yields T* uniquely, it can be shown 
(see Ref. 13) that T* and X^ are continuous at Junction times between 
the velocity-constrained and the unconstrained arcs. Thus , thrust is 
seen to be a continuous function of and t. When C 2 is small, the 
range of X^ for intermediate thrust is narrow; when C 2 ■ 0, the inter- 
mediate thrust i ' singular. (In some cases, not considered in this paper, 
C 2 is negative; then, intermediate thrust is not fuel-optimal.) 
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W« observe thet Bqe. (3) end (IS) for ^ and respectively, and 
the fact that and are «mstant, Is^ly 

ttSO if X^bO onan Interval (22) 

This Is true Irrespective of the olnlalsatlon of R with respect to u. 
From the latter, we obtain, by using ■ 0 and > 0, 

1 W If |u| < axid Xy < 0 

• • 

Ug|8gn|i - %sgnX^ If |p| > Uj, and X^ < 0 

where 

W - •‘gnvX^/(2k2Xv) (23b) 

t ^ And X^ does not vanish on an ir.cerval (denoted by X^ j! 0), 
Che minimization of H gives at once 

u* - -Uj^sgnXjj^ if X^ > 0 and X^ ^ 0 (23c) 

We note that X^ cannot be positive while X^ vanishes on an 
interval because then minimization of H implies u* ■ ±Uq which is 

incompatible with Eq. (22). However, if X^ vanishes on an interval 

and \y crosses zero from negative to positive values, say at t > tj, 
then u switches from u(t) ■ 0, t < tj to u(t^) ■ this is a 
transition from a straight-line flightpath to a curved one. 

We also note that the simultaneous vanishing on an interval of 
X||, and X^ is not fuel-optimal because it implies the vanishing of 
X^ and \y which in turn Implies ” V * nonoptimal 

case of Eq. (18). 

In summary, we see that u, and hence the baidi-engle, are continuous 
in tine when X^ is negative; u is discontinuous when X^ is positive. 
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or, ot tho mamnt ccoomo Mto to boooM poiitivt «Ml had bam 
aarp on tin prevloua interval. In view of Iq. (21), the diaeontl nulty 
of the b ank-ani|la ocpnr a at nlnif tteuat. 

Laotly, ee evalnate p. Whan v(t) s v—>. the thruat la inter- 
awdlate in^r Aaauaptlon 1. Then, ^ Bqa. (S) and (20), T ■ D ^vea 

-(Cj + \v/«)/2Cj " kjV^ + kj/viia ♦ k^«*/v»„ 

This equation gives an expression for and, iqpon differentiation, 
for X^, which we substitute into Bq. (16) to obtain 


n ■ -Xjj cos p 


- + iCjlk^v^ + kj(l + - k^a + 


(24) 

For a velocity-constrained arc to be optimal, n(t) given by Bq. (24) 
ORist be nonnegative. 


IV. Optimality of Straight-Line and Circular Pata Blaaents 

Since the suboptimal algorithms of Refs. 10 and 11 are based on 
piecing together circular arcs and at most one straight-line segment, 
we are interested in the optimality of these path elements. 

We first show that there can be at most one strai^t-line segment 
in a fuel optimal path and, if it occurs, it must do so at the beginning^ 
of the pat^. A subsequent curved path, if any, starts by a switch to 
maximum bank angle Ul * at (continuous) minimum thrus^ ^ “ ^min» 
the bank angle's mignitude and the thrust remain_at £nd 
respectively, a s lo n g as Xy re main s positive. 
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To prove cblo propoeition* we note tbet fU(|it elmip e etralglht-liiie 
eefemt on e eublntervel (t^» t^J ie cherecterlsed Ivy 

♦(t) - • u(t) ■ 0 • X^(t) , Xy(t) < 0 on (tj, tjJ (25) 

He observe thst for a strel^t-Une peth, the point ((> ■ ■ 0) Is 

an equilibrium point* for Bqs. (3) and (IS) for p and with *'(t) 
as a continuous parameter. Hence the straight-line aeipient can be 
entered at t^ and exited at tg by an optimal control only if the 
control is discontinuous at t^ and tj. Therefore, as noted in the dis- 
cussion following Eq. (24), it is necessary that Xy(t) cross sero at 
tj and t 2 according to 



(26) 

Xy(t2) - 0 , ly(t2) > 0 

(27) 


(This causes u(t) to switch from u(t7) “ tUg, to u(t) - 0 on 

(tj, t 2 ] and back to u(t^) ■> tUg,.) Thus, Xy(t^) • 0, 1 ■ 1, 2; by 

Eq. (21) this inqplies T(t^) • Since X^(t) la continuous, we 

have X^(t£> ■ 0. Using X^(tj) « 0 in Eq. (15), and X^(t^) ■ X^(tp ■ 0 

and T(t^) - Ty,^ n in Eq. (17), gives two equations in X^ and Xy which 

yield 


*An equilibrium point of differential equation x - f(x,t) is a 
point x^(t) that satisfies 0 » f[x(t), t]. If f(x,t) is continuous 
in X and t, then a solution x(t) of the differential equation started 
outside (at) the equilibrium point cannot reach it (depart from it) in 
finite tine (see, e.g. , Kef. 14). 
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1 - I, 2 


08) 


X, « CO* 8«>/v(c^> , 


Xy - -(1^ *1II ♦•)/v(4> . 1 - 1, 2 (29) 

i^r* 1 ^ 4 ^ 1* th* fti*l flow r*t« ai; 

^■in • «o *► «x^ ♦ « An <»> 

Using Bq*. (28)-(30) In Bq. (16) giv** 

^v(ti) ■ • 1 ■ 1. 2 (30 

Since by 4seunpclon 2 Is positive, Bq. (31) shows that 

Xy(t|^) > 0, i > 1, 2. This contradicts Bq. (26) but confim Bq. (27). 
Hcoice, a curved flightpath cannot precede, but can follow, a straight- 
line segment. On the curved path, at least initially, X^(t) > 0, t > t 2 . 
Hence, by Eq. (23c), |u(t)| * Ug,, and by Bq. (21), T ■ T^n. This com- 
pletes the proof. 

Next, we ask whether cruise at a constant velocity v(t) s < v,^v 
on a straight-line segment Is fuel optloal. For the purpose of the 
inquiry, we temporarily lift the constraint v(t) < Vg^^. He find that 
a flight at constant velocity on a straight-line segMnt of a fli^tpat^ 
CM be fuel optimal only if the e nti re flightpath is s traight, flown at 
a c onstant o ptimal cruise ve locitjr v^* ^vm by J^he solution of t^ 
equatjUm 


SCjkj^Vc* + CjkjVg* - (Cg 2cjkjkj)Vg^ - 3c,kjVg* - Scjk^ - 0 ; (32) 

Eq. (32) has one and only one real and positive solution. 

Proof . For flight on a straight path at a constsnt velocity v^. the 
thrust is constant, T ■ T^ ■ ^(Vc)t 6y Asstn^tion 1, is intermediate; 
brace, must be constant, , according to Eq. (20). The point 
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(tl* ■■ V ■ ■ 0» Xy ■ Xy ) !• an •qvilibrim point for the respec- 

tive set of differential eqtMtione [Bqs. (3)« (9)» (15) » end (16)], Whence 
folloira the first part of the proposition. To derive Bq. (32), we 
elinlnete Xjj end Xy in Bq. (16) hy using Bq. (17), end we substitute 

the value of Xy ■ Xy from Bq. (20). This yields 

c 

XyVg - {cj ♦ Cj(kjVc* + kj/vg*) + Cj(kjV^* ♦ 

- 2(2c2(kiVc* + kj/vc*) + cJ(kjVc* - kj/vc*)) (33) 

Since X^ is a constant, the right aide of Bq. (33) vanishes, yielding 
Eq. (32). Let the right side of Eq. (33) be denoted by f(v^). Then 

df(v^)/dv^ - -2(Cjk^Vg + ■*" ^OCjkj^/v^®) (34) 

For sufficiently small v^, f(v^) > 0, and from Eq. (34), df/dv^ < 0 
for all positive v^. Hence, f(v^) is monotonically decreasing and has 
one and only one real and positive zero v^. • v^*. This completes the 
proof. 

The significance of the velocity v^* i s that when the enti re 
fligh tpath la strai^£, v^* ^ bo opper botmd on v(t), t c [0, t^) 
in the se nse 

if v(0) < Vg* and V£ < v^* , th«i v(t) < v^* (35) 

and a lower boi^ in th e sense 

^ v(0) > Vg* and V| > Vg* , then v(t) > Vg* (36) 

Proof . Consider Eq. iJS), and stq>pose that the velocity lises to a 

maximum v(t*) > v *; there v(t*) ■ 0, and v(t') < 0 which we want 
to contradict. At t ■ t* we have 

v(t') “ (T - (dD/dv)v(t')J/m - f/m ■ -2Cjiy(t*)/m* (37) 
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Since v(t*) ■ 0, Bq. (33) applies vith v(C') replMlag v^. For v(t*) > 
the left side of Eq. (33) venishee. Since hy Bq. (34) the rigM^ side of 
Bq. (33) is decreasing, then, if v(t') > v^*, we have ly(t*) < 0. It 
follows from Bq. (37) that v(t*) > 0, which is s contradiction, fhis 
proves Bq. (35); Bq. (36) is proved analogously. We note in passing 
that by linearizing the equations for v and around the equilibriun 
point (v.*, ), this point can be shown to be a saddle point, which 

conforms with Eqs. (35) and (36). 

The case of Eq. (36) is of little interest for landing because 
usually V£ < v^*. If v(0) and v£ are on opposite sides of v^*, then 
v(t) can cross v^*. Our numerical experience, described in the next 
section, shows that acts as an upper bound on v(t) also for 

relatively long-range fllghtpaths which are not strictly straight. Con- 
sidering now the constraint v(t) < for the case of Eq. (35), it 

is clear that it will be inactive if V maif ^ v^.* and is likely to be 
active if Vg^^ is much below v,.*. 

The optimal cruise velocity v^.* given by Eq. (32) is the velocity 
that provides minimum fuel consumption per unit distance along a straight 
flightpath. One expects that v^* will be higher if C 2 * 0 in the 
fuel consumption model of Eq. (10), and it will be lower if Cg • 0. 

For our numerical values 

! 349.5 knots if Cp # 0, Cj ^ 0 

359.0 knots if Cq ^ 0, c^ ■ 0 

298.8 knots if Cp ■ 0, c^ “ 0 

For comparison, the minimum-drag velocity, which is the minimum fuel per 

unit time velocity, is 

Vp “ (kj/kj)^^** • 227 knots 
*^mln 
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BxaBlnatlon of the necesaary conditions shows that a circular 


fli^tpa th can ba fuel optimal only if both the baiJt"angta and the 
velo city are at thait respective constraint bound s* The proof is 
straightforward and is omitted. 

V. Computation of Extremals 

The extremals are conQtuted by numerical integration of the state 
and costate equations [Eqs. (l)-(S) and (12)-(16)]» with the controls 
given by Eqs. (19) and (23) • from a given state and with chosen initial 
values of the costace variables as parameters. One can thus obtain 
families of extremals but one cannot meet specified end-conditions. 

In the computations, the state constraint v < v^^ » 250 knots 
was imposed only on some of the long-range flightpaths; however, no 
extremal was extended beyond 350 knots. A state constraint v > v ^^n 
was not explicitly incorporated in the necessary conditions; however, 
extremals where v(t) dropped below 150 knots were rejected because the 
drag model of Eq. (5) would not be valid below that velocity. 

In the following figures, distances are in nautical miles, the 
velocities v are in knots, the thrust T is in thousands of pounds, 
the fuel consumption f to the final point is in pounds, and the time 
in seconds is from the integration's starting time t^ > 0. The portions 
of the fllghtpath with T > 0 are in bold curves. The maximal and 
minimal values of thrust and velocity are among the T's and v's shown. 
The arrows show the airplane's direction of flight. The starting values 
of the costates at tg ■ 0 are shown in the figures by 
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(to o«r eo^tttolims m divitod eh« naailtoiktoB bf c^). 

It is convenient to stert Integretimi bechverds free the knoim and 
fixed final state of Bq. (9). Confuting X^ fro* the condition H *■ 0 
at the integration's starting tinie tg» the extresals are determined by 
the parameters X^. X^(tg), and X^Ctg). Since the problem is time- 
invariant, we set tg » 0. The sign and relative magnitude of both 
Xy and X^(tf) determine the direction of turn; it can be shown that by 
changing the sign of both, the fllghtpath is reflected about the x-axis 
(true if integration is started at y > <|i > 0). Such backward integra- 
tion produced coasting and short-range turning extrenals, as those shown 
in Figs. 1 and 2. 

Figure 1 shows coasting flightpaths, namely, decelerating paths with 
zero thrust throughout. For extremals 1 and 2 the bank angle switches 
betifeen its bounds of ±30° because X^(t) is positive throughout. This 
is typical of most, but not all, coasting extremals: extremal 3 ends 
with a smooth transition to a shallow bank angle, and extremal 4 starts 
with a straight-line segment. Coasting extremals may be significant for 
emergency landing. 

Figure 2 shows (by solid curves) partly thrusting extremal paths 
whose backward computation was arbitrarily terminated at 200 knots. These 
paths represent short-range turning approaches starting at 200 knots, such 
as after an aborted landing. Typically, for a turn through a large angle, 
as in paths 1 and 2, the velocity first drops to achieve a tighter turn. 
This was noted in Ref. 6 for mi’'-imum-time turns, but it is less intuitively 
obvious in the minimum-fuel case because fuel is later expended to 
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Fig. 1 Coasting (zero-thrust) flightpaths. 
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Fig. 2 Short-range turning flightpaths. 
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accelerate. The broken curves* here and In Fig. A* will be connented on 
in the next section. 

Backward integration from the final state did not produce long- 
range fllghtpaths. Such paths require a sustained intermediate thrust 
which is dictated by a narrow range of (from Bq. (20) it follows at 
once that T > « 0 if X^/c^m > -1.0* and by computation* T i 

if X^/c^m < -1.215]* resulting in extreme sensitivity to the choice of 
Xv(tg). Therefore* long-range extremals* and other types of extremals 
with special conditions at intermediate points, were produced by forward 
and backward integration from an appropriate Intermediate state. For 
convenience, integration was started at x(ts) ■ y(tg) ■ <|)(ts) ■ 0, 
tg > 0, with an appropriate v(tg). The resulting flightpath can then 
be shifted and rotated to satisfy the final condition of Bq. (9). 

Integration of long-range extremals without the velocity constraint 
[Eq. (8)1 was started at a velocity below the optimal cruise velocity 
v^* ■ 349.5, with such that T « D. The parameters are then v(tg), 
Xy and X^(tg), tg » 0. The velocity profiles of such paths are shown in 
Fig. 3. Since v(tg) <■ 0, v(tg) is the maximal velocity. Curves 1-3 
correspond to straight-line fllghtpaths. Curves 4 and 5 correspond to 
paths with initial and final turns as indicated. For t < 0, curves 3, 

4, and 5 are indistinguishable. On curve 4, we note the dip in velocity 
at the large final turn to vCt^) • 250 knots. On curve 5 the velocity 
decreases faster than on curve 3 because of the added drag due to maximum 
bank-angle. 

Bxtremals with a velocity-constrained arc were confuted by starting 
on the arc at v(tg) - v^g^ * 250 knots. The costate Xy(tg) is 
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determined by the requirement that thrust equal drag. The remaining 
parameters are end two parameters to control the departures 

from the velocity-constrained arc according to preset times or according 
to the multiplier n(t) >0 computed by Bq. (24). We considered 
extremals with only one constrained arc* of the long-range type. 

Figure 4 shows long-range fllghtpaths* which typically consist of a 
possible initial turn* a long* almost strai^t arc* and a possible final 
turn. Path 1 is a velocity-constrained extremal* path 2 is unconstrained* 
and both were conq>uted from the point indicated by t ■ 0. For convenience 
of presentation, neither path is particularly long. In both cases the 
almost straight arc can be made as long as desired: in the constrained 

case by selecting and X||,(0) sufficiently small* and in the uncon- 
strained case by selecting* in addition, v(0) sufficiently close to the 
optimal cruise velocity v^*. 

Integration from an intermediate state, rather than frcm an end- 
state, is mandatory also for the type of extremal discussed in Sec. IV: 
a straight-line segment followed by a switch to a maximal bank-angle turn. 
This extremal requires the satisfaction of Bq. (31) at the switch-time tj. 
Since (j/(t), u(t), X|p(t), and are zero on the straight-line segment, 
the remaining parameters are v(t„) and X (t.l < 0. For our numerical 
values and range of velocity, these extremals (such as fllghtpath 4 in 
Fig. 1 whose integration started at the point indicated by v ■ 300 knots) 
are all coasting extremals with maximal bank-angle throughout the turn. 

The accuracy of the computed extremals was checked by decreasing the 
integration step and by observing the accuracy of the condition H ■ 0. 
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Fig. 4 Long-range fli^tpaths. 




VI. Optimal and Near-Optimal Solutions 


We compare now the fuel consunq>tion of fUgihtpaths produced by 
extremals with that of near-optimal paths. We thus test the quality of 
the near-optimal paths as well as the global optimality of the extremals. 

The near-optimal algorithm of Ref. 10 generates a flyable state- 
trajnctory between any two end states; it is sufficiently fast to be 
operated on-board and in real time. The algorithm, based on the results 
in Ref. 1, generates the shortest flyable path consisting of circular 
paths Joined by at most one straight-line segment. In Ref. 11, the algo- 
rithm is further refined in the li^t of the results of this study; in 
particular, the curved path is created by a succession of 30” circular 
arcs of varying radii. In view of the results of Sec. Ill, it is clear 
that such a synthesized flightpath cannot be fuel-optimal. However, it 
is evident from Fig. 4 that such paths can closely approximate the 
optimal ones for the most common and important type of path, the long- 
range path. The details of the algorithm as well as numerous comparisons 
are reported in Ref. 11. For example, for 28 long-range paths of about 
20 miles, without a velocity constraint, the worst approximation was 
2.95Z off the minimum fuel consumption, the best 0.5Z, and the average 
1.S2Z (Ref. 11). 

The extremals exhibited in Sec. V satisfy only necessary conditions 
for optimality. Are they optimal? He have no proof, but fuel-optimality 
can be argued for at least the coasting extremals, which are also minimum- 
time coasting extremals. Consider flightpaths 1 and 2 in Fig. 1. Along 
these the magnitude of the bank-angle, and hence the drag, are at all 
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tlMS naxinal. Thar«for«» any other coeating tUghtpaOi with 
x<t{) ■ y(t^) - t(i(t£) - 0, t^ < t£, Mot h«v« v(tp > - 180 knots. 

The fuel-optlnelity of paths 3 and 4 can be supported by alnllar, albeit 
sonatfhat weaker erginents. 

The situation la far leas clear for the thrusting extrenals. These* 
as the integration continued, often started twisting and looping In a 
manner which appears increasingly nonoptimal, as shown by the broken 
curves of Figs. 2 and 4. Such behavior of extrema^) is likely to have 
been observed pr< vlously; there is an allusion to loss of global optl** 
mality in Ref. 7. In general, as an extremal is extended by integration 
from some starting point, a time t^ may be reached beyond which the 
extroaal ceases to be globally optimal; tj) is called a Darboux point 
with respect to the starting point (Ref. 12). 

The existence of Darboux points on the extremals of Figs. 2 and 4 
is doDonstrated in Fig. 5 by comparison of fuel consumption with near- 
optimal paths. Figure 5a shows path 1 of Fig. 4 which, as the forward 
integration is extended, has in the final turn three points with 
v ■ 180 knots. The comparison confirms the nonoptisiality of end-points 2 
and 3; the Darboux point appears to be about midway between points 1 and 2. 
We note that also path 2 of Fig. 4 has two possible end-points with 
V > 180 knots. In this case, however, the extremal path to the second 
end-point is deemed to be optimal. Figure 5b shows the fuel consumption 
to the final point x • y ■ 0 of a turning extremal path and of near- 
optimal paths, for several points along the extremal; three such near- 
optimal paths are shown. Evidently, for points beyond v ■ 320 knots or 
so, the near-optimal paths use less fuel. Altlx>ugh we do not have the 


22 


v-in 



Fig. 5 Comparison of extremal and near-optimal flightpaths. 
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51g. 3 Concluded 



•met location *f tte Darbonx point* na ara confidant that If tha axtrmal 
la tamlnatad at* aay* 2S0 knota* It 5s optiml. 

Tliua* althouih ona cannot naa tha naar-optlnal patha to prova optl- 
Mlltp* ona can gat a roiq^ Idaa of tha location of tha Daihonx point* 
particularly for tha long^ai^ patha* utera tte naar-opttoal an>roxlma* 
tlon la vary good. Of couraa* portlona of nonoptinal axtramla nay ba 
optimal; for axamvlo* tha coaatlng portlona of tha axtramala In Pig. 2 
are llkaly to be optimal. 

VII. SuuMry md Concluding Ranarke 

The characteristics of minimum fuel horlsontal flightpaths in the 
terminal area were investigated analytically and computationally. Analysis 
of the necessary conditions showed the following. 

1) Thrust is continuous, but the baidcangle may be switching for 
certain values of the costates Xy and X^ (see Eq. (23c)]. 

2) A straight-line segment may be fuel-optimal only if it is at the 
beginning of the flightpath; a subsequent curved path* if any* must start 
by a switch to maximum bank angle while the thrust is at minimum value. 

3) The optimal cruise velocity* given by the solution of Eq. (32)* 
acts as an upper bound on the velocity for straight (or almost straight) 
minimum-fuel flightpaths. 

4) A circular flightpath may be fuel optimal only if both the bank 
angle and the velocity are at their respective constraint bounds. 

The coiq^tation of extremals produced many representative minimiBB- 
fuel flightpaths that can be categorlaad as Im^ranga paths* short-range 
turning paths* and coasting (aero-thrust) patiM. US found that: 
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1) Extreme sensitivity to clioice of for long-range paths could 
be overcome by starting totegratioa at an appropriate intermediate state* 
rather than an end-state. 

2) Long-range paths with a large initial turn (over 100*) start by 
deceleration followed by an acceleration in the remainder of the turn. 

3) Long-range paths with a final turn up to 140* end with zero 
thrust and turn with maximum baidc-angle magnitude. However* if the final 
turn is large* and in particular if in addition the final velocity vCt^) 
is higher than 180 knots* the turn is executed by decelerating below 
v(tf) and a final acceleration at maximum bank-angle magnitude. It is 
shown in Ref. 11 how these findings made possible the refinement of an 
existing on-line algorithm of Ref. 10 to the point where the fuel consump- 
tion of long-range near-optimal paths is well within 1-31 of that of 
optimal paths. 

The near-optimal algorithm was very helpful also in alleviating the 
problem of finding the Darboux points. We found that: 

1) Turning extremals that require thrust toward the end of the path* 
produced by backward integration from v(tf) « 180 knots* inevitably 
became nonoptlmal at some point beyond v ■ 250 knots. 

2) Turning but coasting extremals* on the other hand* appear to be 
optimal no matter how long they were extended (backwards from 

v(t£) B 180 knots) . 

3} Optimality of long-range extremals may be lost if integration is 
extended so that the Initial or final turns are much larger than 180*. 

Evaluation of the optimality of extremals by a near-optimal algorithm 
proved to be a practical solution to the Darboux point question for 
flightpath problems. The general problon* however* is unresolved, for 
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there is no test for Darboux points. Could a test (e.g. , Ref. 15) for 
conjugate points (beyond which an extremal ceases to be locally optimal) 
be helpful? The answer appears to be negative; certainly so for a special 
two<-dlmenslonal case of the present problem, in which x(t^) Is free 
and v(t) B constant, which we examined In detail In Ref. 16. The ques- 
tion of global optimality, highlighted here by computation of extremals. 

Is inherent (though perhaps less visible) In other optimization techniques. 
The problem of Darboux points remains a challenge for future research. 

The work reported here is of course but one element in the develop- 
ment of a practical, fuel-efficient, and safe system for terminal opera- 
tion. In particular, an extension to include the third dimension, 
altitude, is to be studied. 

Appendix A: Equations of Motion 

The lateral, longitudinal, and vertical force equations are, 
respectively, 

mv^ B -(L + T sin o)sin ^ (Al) 

mv = T cos a - D (A2) 

mg B (L -I- T sin a)cos ^ (A3) 

For small angle of attack a, Eqs. (3) and (4) result, where u » tan 

Lift and drag are given by L > Ct aqS, D b Cq qS eLa, where c is 

"o o 

the efficiency factor and q « (l/2)pv^ Is the dynamic pressure, S the 

wing area, and the coefficients Ct and Cn are assumed to be independent 

a o 

of the velocity. Now, 

eLo b eL^/(CL qS) - e[mg/cos 4 - T sin a]*/(CL qS) 
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Neglecting T sin a, D is of Che form of Eq. (5), with • (l/2)(^ pS» 

o 

k. * 2em^g^/(CT pS). Using Che values €{> • 0.01S» C|. • 0.08/deg, 

* "O o o 

e - 0.004/deg, (l/2)p - 295~' lb/ft*/knot*, S - 1560 ft*, mg - 150,000 lb, 
gives Che values of and k 2 in Sec. II. 


28 


Acknotfledgnents 


The authors wish to thank l>r. Heinz Brsberger for fruitful discussions 
and perceptive couDients. This research was supported, for the first 
author, by a National Research Council Assoclatesbip at bines Research 
Center, NASA. 


References 

^Erzberger, H. and Lee, H. Q. , "Optimum Horizontal Guidance Techniques 
for Aircraft," Journal of Aircraft , Vol. 8, No. 2, Feb. 1971, pp. 95-101. 

^Pecsvaradi, T. , "Optimal Horizontal Guidance for Aircraft on the 
Terminal Area," IEEE Transactions on Automatic Control , Vol. AC-17, 

No. 6. Dec. 1972, pp. 763-772. 

^Kishi, F. H. and Pfeffer, I., "Approach Guidance to Circular Flight 
Paths," Journal of Aircraft , Vol. 8, No. 2, Feb. 1971, pp. 89-95. 

''Hedrick, J. K. and Bryson, A. E., Jr., "Minimum Time Turns for a 
Supersonic Airplane at Constant Altitude," Journal of Aircraft , Vol. 8, 

No. 3, Mar. 1971, pp. 182-187. 

^Hoffman, W. C. and Bryson, A. E. , Jr. , "Minimum Time Turns to a 
Specified Track," ASI-TR-71-4, Aerospace Systems, Inc., Sept. 1971. 
Summarized in: Hoffman, W. C. , Zvara, J., and Bryson, A. E., Jr., 

"Optimum Turns to a Specified Track for a Supersonic Aircraft," Preprints 
of Joint Automatic Control Conference, 1972, pp. 955-956. 

^Bryson, A. E. , Jr. and Parson, M. G. , "Constant Altitude Minimum 
Time Turns to a Line and to a Point for a Supersonic Aircraft with a 
Constraint on Maximum Velocity," SUDAAR No. 437, Stanford Ihilverslty, 
Stanford, Calif., Nov. 1971. 


29 



^Hoffman. W. C. and Bryson, A. B. , Jr., "Minimum Time Maneuvers to 
Specified Terminal Conditions," ASI-TR-73-12, Aerospace Systons, Inc., 

Jan. 1973, Sect. 3.2. 

^Bryson, A. E. , Jr. «id Lele, M. M. , "Minimum Fuel Lateral Turns at 
Constant Altitude," AlAA Journal , Vol. 7, No. 3, Mar. 1969, pp. 559-560. 

^Vlnh, N. X., "Mlniffltim Fuel Rocket Maneuvers in Horizontal Fli^t," 
AIAA Journal , Vol. 7, No. 2, Feb. 1973, pp. 165-169. 

^®Erzberger, H. and McLean, J. D., "Fuel Conservative Guidance System 
for Powered Lift Aircraft," Journal of Guidance and Control , Vol. 4, 

Mo. 3, May/June 1981, pp. 253-261. 

^ ^Neuman, F. , "An Algorithm for On-line Construction of Fuel Conserva- 
tive Capture Trajectories," in print. 

^^Mereau, P. M. and Powers, W. F. , "The Darboux Point," Journal of 
Optimization Theory and Applications , Vol. 17, Nos. 5/6, 1975, pp. 545-559. 

^’Jacobson, D. H. , Lele, M. M. , and Speyer, J. L. , "New Necessary Con- 
ditions of Optimality for Control Problems with State-Variable Inequality 
Constraints," Journal of Mathematical Analysis and Applications , Vol. 35, 
1971, pp. 255-284. 

^“Stem, T. E., Theory of Non-Linear Networks and Systems. An 
Introd\iction , Addison-Wesley Ptib. Co., Inc., Reading, Mass., 1965. 

^®Moyer, H. G. , "Optimal Control Problems that Test for Envelope 
ontacts," Journal of Optimization Theory and Applications , Vol. 6, No. 4, 
1970, pp. 387-398. 

^®Kreindler, E. and Neuman, F. , "Global Optimality of Extremals: An 

Exmi^le," Journal of Optimization Theory and Applications , Vol. 37, No. 4, 
Aug. 1982 (to appear). Also NASA TM-8 1240, 1980. 


30 



PART 2t THE CASE OF SINGULAR THRUST 


Elleser Rretndler 


Introduction 

In Ref. 1, the fuel flow rate model Incliides a amall term quadratic 
in thrust. When this term is neglected, as is usual, the intermediate 
thrust is singular. In this note the singular thrust is derived, and 
the Generalized Legendre-Clebsch condition Is examined for various ranges 
of bank-angle and velocity. The results are summarized in the last sec- 
tion. The literature on flightpath optimization in the horizontal plane 
is reviewed in Ref. 1. 


Problem Formulation 

The point-mass equations of motion in the horizontal plane are 


V cos A 

(U) 

V sin A 

(lb) 

-gu/v 

(U) 

(T - D)/m 

(Id) 


Here x and y are the coordinates of the horizontal plane, <|i the heading 
angle measured counterclockwise from the x-axis, v the velocity, g 
the gravitational constant, and m the mass. The control variables are 
the thrust T and u, where u “ tan ^ and A is the bank-utgle, positive 
with right wing down. The drag D is given by 

D ■ kjV* + kj(l + u*)/v* (2) 
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where and k 2 are constants. These equations are derived by assuming 
aero wind velocity, constant mass, coordinated turns, and a small an gio of 
attack which is automatically adjusted to maintain horizontal flight (Ref. 1) . 
The controls are constrained by 

S <Hi, (3) 

and the fuel consumption to be minimized is given by 

J - f ^ ICo + CiT(t)Jdt (5) 

•o 

where c, and c^ are constants and t£ is free. In Ref. 1. the fuel 
flow rate model includes a quadratic term C 2 T^(t) with a small C 2 . 

Application of the Minimum Principle 


The Hamiltonian is given by 

H ■ Cq + CjT + XjjV cos ill + XyV sin iji - X^gu/v + Xy(T - k^v* - kj(l + u^)/v^]/m 


The costate variables are given by 

(6) 

Xjj ■ -Hjj ■ 0 =» X^ ■ constant 

(7a) 

Xy ■ -Hy ■ 0 •» Xy ■ constant 

(7b) 

X^ ■ » v(Xjj sin i|» - Xy cos i|») 

(7c) 


Xy - -Hy ■ -Xjj cos i|» - Xy sin (|> - X^gu/v 

+ 2X^lkjV - kjCl + u*)/v»l/m (7d) 
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where ■ BH/Sx, etc. Since is not prescribed end ■ 0. we 
have the necessary condition 

H - 0 , all t e [0, tj J (8) 


The Hamiltonian is linear in T, with the term multiplying T given 


by 


Hj ■ Cl + X^/m 


(9) 


Minimization of H with respect to T gives 


T* 


'max if 

f«in if » 


the thrust can be intermediate only In the singular case 
Hj, « 0 on a subinterval of [0, t£] 


( 10 ) 


In this note, we are concerned only with the case (10); therefore. 


Xy - -Cjm 


Since in this note X^ is negative, minimization of H with 
respect to u yields 


u* 


gvX^/(2c^k2) if |gvA^/(2c^k2)| < u„ 

u„8gnX^ if lgvX^/(2cik2)| > u„ 


(lla) 

(llb) 


The singular case (10) implies the vanishing of all time derivatives of 
Hj. Let the first time derivative of in which T appears explicitly 
be the 2qth (it is always even); q is the order of the singular arc. 

The Generalized Legendre-Clebsch necessary condition requires (Ref. 2) that 



The Singular Thrust 


Consider first the case |u| < Ug. Using Bqs. (16) and (18a) in 
^ « 0 ■ Xy. we compute to be 

Hi ■ “8(^x c®® ♦)«/»» - 2cj(kjv'* + 31 t 2)(T - D)/m*v^ 

First, we note that 

(Hj)^ - -2ci(kiv'* + 3k2)/m*v** < 0 

so that the Generalized Legendre-Clebsch necessary condition of Eq. (12) is 

•« 

satisfied. Second, setting Hy > 0 yields the intermediate, singular 
thrust: 

mgy(Xj^ sin ■* Xy cos t|»)u 

T - D (13) 

2c^(kjV* + Ik^/v^) 

Adding the condition H-0 to we obtain 

X^ - ±2c/kj(3kj - kjV** + CjV*/Cj) /gv 
and hence, using Eq. (11a), 

u - it^Okj - kjV^ + c,y 2 /Cj)/k 2 (14) 

Using Eq. (14), we find that the velocities for which |u| < Ug, are 
outside the range [v^, V 2 ], where and V 2 are the real and positive 
roots of 

kjv** - CgvVc^ + k2(u„,* - 3) - 0 (15) 

(Of course, real and positive v^ or V 2 , or both v^ and V 2 , may not exist; 
in particular, if U|g^ > 3 and (cg/c^)^ < 4k2k2(Ujg^ ** 3) , then |u| < Ujg 
for all V.) 
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We observe that for a straight-line flight* aettlng u ■ 0 In 
Bq. (13) gives T > D. which Implies v(t) ■ - constant. Setting 

u « 0 In Eq. (14)* we obtain the optimal cruise velocity v^* 

(Vc*)* - Co/(2kiCi) + I(co*/(2kiCi)* + akj/kxl'^* (16) 

We now consider the case |u| ■ which occurs when the velocity 
enters the range (v^* V 2 ]. Proceeding as In the previous case for 
|u| < ti|Q* we obtain 

(ftr)^ - -3(CjD - Co)/(mv)* , (17) 

2mgv(X sin - X cos i|))u-8gnX^ 

T - D i 2- , (18) 

3(CjD - Co) 

and 

■ Ci[3k2(l + %*) - kiv" + CoV^/cil/(2gvUj,) 

We observe from Eq. (17) that the Generalized Legendre-Clebsch 
condition of Eq. (12) is satisfied as long as 

D - kjV* + k2(l + Ubj^)/v^ > Cj/Cj 
or 

kjv' - c,v*/Cj + kj(l + u„*) > 0 (19) 

The left side of Eq. (19) has no real zeros, and hence Eq. (19) holds if 

CoV(2cj)* < kjkjd + u„*) (20) 

If Eq. (20) is violated (as when Cq Is increased to create a cond>ined 
time and fuel cost functional)* then there exists a range of velocities 
[v|* Vjl* where v^ and v^ ate the real and positive zeros of the left 
side of Eq. (19)* for which the singular thrust of Eq. (18) is nonoptlmal. 
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We observe that [v^, Vj] c: [v^* V 2 lt and therefore the Generalized 
Legendre-Clebsch condition, is satisfied whenever v is such that 

l«l < %• 

We note that letting u ■ u^, T in Eq. (13) is different from that 
in Eq. (18) (the author is indebted to Dr. H. Erzberger for this observa- 
tion). This shows that the intermediate thrust is discontinuous at a 
time ti, say, when the bank angle saturates. This is because 
contains u which is discontinuous at t^ and so must be T(t 2 ) in 
order to satisfy = 0. 

Summary of Results 

For Intermediate bank-angle, the intermediate singular thrust is 
given by Eq. (13) and it satisfies the Generalized Legendre-Clebsch 
condition. The bank-angle is given (except .for sign) by Eq. (14) and 
the velocity is outside the range [v^, v^], with v^ and v^ being 
the real and positive roots of Eq. (IS). On a straight-line flight, 
the singular thrust is constant and the constant velocity is the optimal 
cruise velocity v^* given by Eq. (16). 

When the velocity enters the range [v^, Vj] the bank-angle 
saturates, |u| - Ug,, and the intermediate singular thrust, now given by 
Eq. (18), undergoes a Jump. The singular thrust satisfies the Generalized 
Legendre-Clebsch condition if Eq. (20) bolds; if Eq. (20) does not hold, 
there exists a range of velocities [v[, v^l c [v^, Vj], where v[ and V 2 
are the zeros of the left side of Eq. (19), such that intermediate thrust 
is nonoptimal. 
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